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ABSTRACT 

High-redshift gamma-ray bursts (GRBs) offer an extraordinary opportunity to study aspects of the early 
Universe, including the cosmic star formation rate (SFR). Motivated by the two recent highest-z GRBs, GRB 
080913atz ~ 6.7andGRB 090423 at z ~ 8.1, and more than four years of ^vvZ/f observations, we first confirm 
that the GRB rate does not trace the SFR in an unbiased way. Correcting for this, we find that the implied SFR 
to beyond z = 8 is consistent with LBG-based measurements after accounting for unseen galaxies at the faint 
end of the UV luminosity function. We show that this provides support for the integrated star formation in the 
range 6 < z < 8 to have been alone sufficient to reionize the Universe. 

Subject headings: gamma rays: bursts — galaxies: evolution — stars: formation 



1. INTRODUCTION 

The connection b etween gamma-ray bursts^ and core- 
collapse supernovae jStanek et al.l 120031: IHiorth et alj 120031) 
tells us that, in observing a GRB, we are witnessing the 
death of a massive, short-lived star. The intense brightness of 
gamma-ray bursts gives hope that, starting from this principle, 
we can probe t he history of s tar formation to very early times 
jjotani' 1997; Wijers et al. TW, 'Lamb & ReiclMd HQQS 
[Porciani&Madau 2001; Bromm & Loeb 2002), potentially 
to higher redshifts than with galaxies alone. First, we must be 
able to observe the GRBs and obtain redshifts for a sufficient 
number of events. Second, we need to understand how to 
calibrate the GRB rate to the star formation rate (SFR). Swiffi 
jGehrels et al. 2004) has pushed the former greatly ahead, and 
allowed studies of the latter 

Our goals are to use the large set of Swift gamma- 
ray bursts with known redshifts (see Fig. [U accumulated 
over the last > four years to examine the above two 
points in greater detail. With improved statistics, we con- 
firm the finding that gamma-ray bursts are not unbiased 
trac ers of the SFR, as in [K istler et al. (2008) (also see, 
e-g.JDaigne et al.'2006VLe & Derm en2007uYuksel & Kistleri 
12007; Salvaterra & Chincarini 2007), and comment on its sus- 
pected origins. This does not, however, prevent a study of the 
amount of high-z star formation; it in fact allows for a more 
proper estimation. 

Several recent high-z bursts, most notably GRB 080913 at 
z ~ 6.7 jG reiner et al. 2009) and GRB 090423 at z ~ 8.1 
(iSalvaterra et al..2009:,Tanvir et al.l2009l). also allow u s to ex- 
tend the SFR determinations from lYiiksel et al.l (I2008h (which 
went to z 6) to even higher redshifts. Here, direct SFR 
measurements are quite challenging, particularly at the faint 
end of the galaxy luminosity function, where GRBs may be 
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ideal tracers. Even with only several events, we determine 
that the SFR declines only slowly from z ~ 4 to z > 8. This 
may confirm that a substantial amount of star formation oc- 
curs within faint galaxies, in agreement with extrapolations 
of Lyman Break Galaxy (LBG) measurements, and suggests 
that stars may be responsible for cosmic reionization. 

2. GRB SAMPLE 

It is easy to understand, with the combination of uncertain 
extinction corrections, cosmic variance, and selection biases, 
why measurements of the SFR at high redshifts are difficult 
endeavors. Principal among these is that flux-limited surveys 
observe the bright end of the galaxy luminosity function (LF) 
and must correct for the faint end, where much of the star for- 
mation may be occurring. The use of gamma-ray bursts as a 
star formation measure will have its own systematic effects; 
however the opportunity presented to examine very-high red- 
shifts, and possibly unseen faint galaxies, is great, with no 
known backgrounds for a bona fide GRB. 

To calculate the expected redshift distribution of GRBs, we 
combine the comoving GRB rate, nGRB(-z) = ^{z) x p*(z), 
where p^, (z) is the SFR density and £(z) accounts for the frac- 
tion of stars resulting in GRBs, with the ability to observe 
the GRB and obtain a redshift (0 < F{z) < 1), the frac- 
tion of GRBs unobserv able due to bea ming ((/beam); e.g., 
IRacusin et al] 120091: ICenko et all l2009h . and the comovmg 




Fig. 1. — The luminosity-redshift distribution of 119 Swift GRBs, as we 
determine from the (updated) Butler et al. (2007) catalog. Squares represent 
the 63 GRBs used in Yiiksel et al. 1 2008), with 56 found subsequently: before 
(grey circles) and after (red circles) the start of Fermi. Three Fermi-hAI GeV 
bursts (triangles) are shown (but not used in our analysis). The shaded region 
approximates an effective threshold for detection. Demarcated are the GRB 
subsamples used to estimate the SFR. 
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volume per unit redshift"^ as 



dN 
dz 



F{z) 



£{z) p^,{z) dV/dz 

(/beam) 1 + 2 



(1) 



F{z) can be kept constant by considering only bursts with 
luminosities sufficient to be viewed within an entire redshift 
range ( Kistler et al. 2008). We then write £{z) =£oil + z)°', 
with £q a (unknown) constant that conv erts p*(z) to a GRB 
rate (in a given GRB luminosity range ). IKistler et all (l2008h 
found that a = Q (directly tracing the SFR) was inconsistent 
with the data at the 95% level, which favored a ~ 1.5. As 
shown in Fig.lT] many more GRBs have since been detected, 
warranting a reexamination of this result. 

Our sample includes the 63 GRBs used in IKistler et al.l 
(^008) (up to 15 May 2007), supplemented by 56 subsequent 
Swift events with redshifts. We calculate the intrinsic (av- 
eraged) GRB luminosity, Liso = £'iso/[79o/(l + z)], from 
the rest-frame isotropic equivalent (uncorrected for beaming) 

1 — IC* keV energy release (E'iso) and Tgo, the time interval 
containing 90% of the prompt emission, as given in the cat- 
alog'° of Butler"etal] (|2007|) ." The results for GRBs with 
Tgo > 2 s are shown in Fig. [T| 

For this test, we use the cuts defined in lKistler et al.l (120081) : 
GRBs in the range z = - 4 with Liso > 10^^ erg s"^ 
This removes many low-z, low-Ljso bursts that could not have 
been se en at higher- 2;,_leaving us with 66. The SFR fit from 
[Hopkins & Beacom (2(306) in this range is used as a baseline 
for comparison, as shown by the dashed line in Fig. |2] A 
Kolmogorov-Smirnov test confirms that the GRB rate is in- 
compatible with the expectations from the SFR, now at the 
~ 99% level (possibly higher due to likely missing bursts at 

2 ;$ 4) with the present greater statistics, requiring an en- 
hanced evolution relative to the SFR. Even if we exclude the 
range z = 1.5 — 2, where a larger fraction of redshifts might 
be missed (Bloom 2003), the value remains at 98%. Possi- 
ble ori gins of this trend are discussed in detail in Kistl er et al.l 
( 20081). including an ove rall decrease in cosmic metallicity 
dLanger & Normanll2006l) . 

Our present result suggests a slightly lower value of a}^ 
Irrespective of the origins of this bias, it must be accounted 
for in z = 1 — 4 to properly relate the GRB rate to the SFR. 
In Fig. 12] we show a shaded band bounded above by a model 
using a = 0.6 and below by a = 1.8, which can be excluded 
at > 84%. To be conservative, we will assume this evolution 
continues to higher z, considering a = 1.2 throughout. 

3. THE HIGH-Z STAR FORMATION RATE 

We briefly review the framework laid out in lYiiksel et al.l 
( 120081) for calibrating a GRB-based estimate of p,. This is 
based on using GRB and SFR measurements in z = 1 — 4 as 
benchmarks for comparison with bursts of similar luminosity 
in a higher-z range. Using only GRBs that could have been 
detected from anywhere within the volume allows for needed 
empirical calibration, since neither the conversion from GRB 
rate to SFR nor the GRB luminosity function are known a 



' dV/dz = 4Tr(c/Ho)dl{z)/y/(J 
comoving distance, Q 



z)^ flm + t^A. where dc is the 
0.3, f2A = 0.7, and Hq = 70 kin/s/Mpc. 
Updated at http://astro.berkeley.edu/~nat/swift 
' ' Note that these values are estimated based upon measurements made by 
Swift (in the 15 — 150 keV energy band), as explained in Butler et al. (2003). 

This may be due in part to the rate of GRB observations at higher red- 
shifts decreasing noticeably in the period following the cutoff date for our 
initial GRB sample (for reasons unknown). Fortunately, as can be seen in 
Fig- El high-z GRBs detections have since increased (denoted as the period 
after the start of Fermi operations). 
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Fig. 2. — The cumulative distribution of the 66 Swift GRBs with Liso > 
10^^ er g s~^ in 2 = — 4 (solid), as compared to the expectations from the 
SFH of lHopkins & BeacomI I20Q6I) alone (dashed) and additional evolution 
of the form (1 + z)^'^ (dotted). Outside of the shaded region (bounded by 
models with a = 0.6 and 1.8) corresponds to an exclusion of > 84%. 

priori. Part of the challenge is in determining the detection 
threshold versus z, since Swift was designed to maximize 
GRB detection, though not necessarily in a way well-defined 
for our purpose (Band 200^. We show in Fig.[T|an estimated 
threshold (cx df, see Kistl er et al.11200 8') based on the GRB lu- 
minosities, which acts as a guide to make cuts that maximize 
statistics and minimize potential "missing" bursts. 

The cuts and resulting subsamples used for the SFR analy- 
sis are shown in Fig. [l]'^ Fig. [3] shows these in comparison 
to the distribution of Liso values for bursts in z = 1 — 4. 
Bursts within each set will be compared to GRBs within the 
range z = 1 — 4 above the given luminosity cut. We em- 
phasize in advance that, although the final bin contains only 
GRB 090423 (at z ~ 8.1), even this single event is signif- 
ican t, as it would b e quite unlikely if was too low (see 
also lSalvaterra et al] |2009). With this GRB we are entering a 
regime where the age of the star is becoming non-negligible 
compared to the age of the Universe, so we extend this bin to 
z = 8.5 to cover a plausible range in progenitor lifetime. 

The "expected" number of GRBs in z = 1 — 4 is 



M^Z^^t^ I dzF{z)£{z) 



.An 

An 

-A I dzp^{z){l + z) 



p^,{z) dV/dz 

(/beam) 1 + Z 

^ dV/dz 
1 + z ' 



(2) 



in which A — At A57£o -fo/47r(/beam) depends on the 
observing time (At), sky coverage (Ail), and luminosity 
range of GRBs under examination. From the average SFR, 
{p* )zi-z2' same can be performed for the other ranges as 



A 



1 + z 



Our interest is in finding {p*) ^-^^z^ dividing out A (using 
Eq.O. Taking the measured GRB counts, Mf^t^^^ to be rep- 
resentative of the expectations, A/'^^^i'^^' ^'^'^ 



Mitzjldz'^p.[z){x + zr 



{p. 



Ml 



obs 



(4) 



Note that the decrease of {dV/dz)/{l + z) at z > 1.5 (as 

Note that we exclude GRB 0601 16 (not shown), with a possible photo- 
metric redshift of z = 6.6 (as listed in , Butler et al..2007. .) 



The Highest-z GRBs and the SFR 
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shown in Fig. 1 of IKistler et al.l 120081) gives progressively 
more weight to each observed higher-z GRB. 

We show our new determinations of the high-z SFR in 
Fig.|4](assuming a Salpeter 1955 IMF). Error bars correspond 
to 68% Poisson confidence intervals for the binned events 
([Gehrels 1986). We also show as a shaded band the values 
obtained for different assumptions of a, bounded above by 
a = 0.6 and below by a = 1.8, which yields an uncertainty 
smaller than the statistics in the last bins. Variations due to 
changing the Li^o cutoff can be determined from Fig.|3] which 
will typically be less than the statistical uncertainties. We 
have been generally conservative and have also verified that 
using another luminosity estimator, the peak isotropic equiva- 
lent luminosity, yields similar results. Other effects, including 
the selection of z-ranges and the inc lusion/exclu sion of partic- 
ular bursts, are discussed in Yiikse l et al.l (120081) . We mention 
only that none of these affect the basic point that the SFR must 
be large enough to produce the observed GRB counts. 

Depending upon the source of the evolution, our bias cor- 
rection may be unduly underestimating by a factor of 
a few at higher z. The most likely astrophysical explana- 
tion is due to metallicity. GRBs are fou nd to favo r metal- 
poor (Stanek et al. 2006), sub-L* galaxies (iFynbo et ani2003l : 
iLe Floc'h et al. 2003 ; Fruchter et al.,.2006l ), so having a larger 
fraction of the SFR within such hosts would result in a higher 
GRB rate. This could be the case with a steepening faint-end 
slope of the galaxy LF, so that more of p» arises from be- 
low L* (L^ as defined at z). This has been observed between 
z = 0andzw2-3 (see Fig. 7 of Reddy & Steidel 2009). 

While our result at z = 4 — 5 is in basic agreement with ear- 
lier measurements, at the highest-z ranges, LBG studies probe 
only the brightest galaxies and must estimate the faint end of 
the UV LF based on limited data. Our results diverge from 
these if corrections for unseen galaxies are not made. For 
example, we focus upon the measurements in lBouwens et al.l 
((2007, 2008), which are reported (lower triangles in Fig. |4]l 
for an integration down to 0.2 i*^3 (with their adopted dust 
corrections). Fully integrating their UV LFs (which can be re- 
garded as giving a maximum), with faint-end slopes of -1.73, 
-1.66, -1.74, -1.74 for (z) = 3.8, 5.0, 5.9, 7.3, respectively, 
yields the upper set of triangles. 

Within the uncertainties, even the highest redshift fully- 
integrated point now agrees reasonably well with our results, 
and the preference of GRBs for faint galaxies (although the 
exact relation between GRB hosts and star forming galax- 
ies as a whole remains to be determined). We note that the 
IBouwens et al.l ( 1200 8) LF slope at (z) = 7.3 was taken to be 
the same as at (z) = 5.9. If the slope is actually steeper (e.g., 
lYan & Windhorst! I2004i) . then these measurements could be 
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Fig. 3. — The cumulative distribution of Ligo for GRBs in the range z 
1 — 4. Shown are the cutoffs used for our GRB subsamples (as in Fig.[TJ. 
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Fig. 4 . — The cosmic star formation history. Shown are the data com- 
piled in IHopJdns & Beaconi l (2006) (light circles) and contributions from 
Lya Emitters (LAE) |6ta et al. 2008). Recent LBG data is shown for 
two UV LF integrat ions: down to 0.2 L*_^ (down triangles; as given in 
IBouwens et aU2 008!) and complete (up triangles). Our (bias-con'ected) Swift 
gamma-ray burst inferred rates are diamonds, with the shaded band showing 
the range of values resulting from varying the evolutionary parame ter be- 
tween a = 0.6 — 1.8. Also shown is the critical p, from Madau et al.] < 19991) 
for C//csc = 40, 30, 20 (dashed lines, top to bottom). 

higher, although it is difficult to draw definite conclusions, 
due to the limited statistics and uncertainties in dust correc- 
tions (e.g.. Chary et al. 2005). 

4. IMPLICATIONS FOR REIONIZATION 

Transmission in the Gunn-Peterson troughs of high-redshift 
quasars implies that reionization must have been accom- 
plished before z = 6 (Fanetal. 2 006|). AGN seem to 
be insufficient for th is purpose (iSrbinovskv & Wvithe 2007t 
iHopkins et al.ll20()8h . leaving stars as the leading candidate. 
To address t he ability of an obse rved population to reionize 
the Universe, IMadau et al.l (11999 ^ provided an estimate for the 
required SFR to balance recombination, pc, which depends 
upon the fraction of photons that escape their galax y (/f,sr,) 
and th e dumpiness of the IGM (C), updated in Paw lik et al.l 
(120091) as 



Pc{z) 



0.027 MqC/U 
Mpc^ yr 30 



1 + z 



7 



0.0465 



(5) 



For comparison with our empirical SFR, we show in Fig. |4] 
curves of pc as a function of z for C / /osc = 40, 30, and 20. We 
find that our SFR estimates can exceed the required to keep 
the Universe ionized at redshifts as high as z > 8. However, 
this criterion refers to an instantaneous equilibrium, and so 
does not address the requirement that the integrated number 
o f ionizations exceed t he n umber of hydrogen atoms. 

IBolton & Haehnell (12007) estimated the ionizing emissiv- 
ity at z ~ 5 from the Lya forest, and looked at simple mod- 
els of the reionization history under different assumptions for 
the evolution of the ionizing photon emissivity at z > 5. 
While their estimate of emissivity from the measured ion- 
ization rate is sensitive to the calculation of mean-free path, 
IBolton & Haehnelt ( 2007) reached the strong conclusion that 
reionization must have been an extended "photon starved" 
process, and that an emissivity which was constant towards 
higher z would have been insufficient to reionize the Universe 
by z ^ 6. This implies that the ionizing emissivity must have 
been higher prior to the end of reionization than just after its 
conclusion. The origin of this higher emissivity could lie in 
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an increase in one or all of the SFR, the escape fraction, or 
the fraction of massive stars in the IMF. Inspection of Fig. |4] 
suggests that the SFR is as large at or could even be higher 
at z ~ 8 than at z ^ 6, implying that the ionizing photon 
emissivity may not be falling towards redshifts greater than 
z ^ 6. Since both ionizing photons and GRBs are produced 
by massive stars, estimates of the ionizing photon emissivity 
from the GRB rate should be fairly robust against uncertain- 
ties in the IMF at t he high-mass end (or at low masses, e.g., 
IWilkins et alJl2008h . 

We are therefore motivated to ask whether we have ob- 
served enough star formation at z > 6 to reionize the uni- 
verse. To answer this question, we make a simple estimate, 
calculating the number of ionizing photons produced prior to 
z 6 given the observed SFR. For a Salpeter (1955) IMF 
and a metallicity of 1/20 Solar, ~ 4600 ioni zing photons are 
produ ced per baryon incorporated into st ars (iBarkana & Loebl 
120011) (further details are given in Wvit he et al.ll2009l) . Taking 
this value, together with a constant SFR for a time interval 
At, we find the number of photons produced per hydrogen in 
the IGM as 

A/- ^ 4 (l^] ( P- ^ ( 

^ VO.l/ VO.lMoMpc-^yr-V V400Myr 

(6) 

Given our SFR, this impUes that A/"-,, ^ St^ 5(/csc/0.1). 

In order to reionize the Universe, more than one ionizing 
photon per baryon is required to compensate for recombi- 
nations in the ionized IGM. Wvithe & Cen (2007) modeled 
the reionization history including evolution of the clumping 
factor with the restriction that reionization end at z ~ 6. 
These models yielded A/^ 4 at z = 6 under a range of as- 
sumptions for the redshift range and efficiency of Population- 
Ill star formation. This is within our estimated range, pro- 
vided that the escape fraction of ionizing radiation is of order 
10%. This value of escape fraction lies in the range found by 
ISrbinovsky & Wyithd (120081) . who combined semi-analytic 
models of reionization and the galaxy LF with simulations 
of the transmission in the high-z Ly-a forest. While not the 
final word, our results may thus indicate that stars produced 
enough ionizing photons in the range 6 < z < 8 to reionize 



the Universe. 



5. DISCUSSION & CONCLUSIONS 



With the discovery of the first astrophysical source at z > 8, 
Swift has enabled GRBs to realize their potential as beacons 
from the distant past, both into the epoch of reionization and 
in adequate numbers at lower redshifts to allow for sensible 
use of the most remarkable events. Using this wealth of data, 
we have estimated the star formation rate at the earliest times 
yet possible, showing that the star formation rate can remain 
high up to at least z 8. From this, it is plausible that the 
level of star formation was sufficient to reionize the Universe. 

The agreement with direct observations, corrected for 
galaxies below detection thresholds, suggests that our GRB- 
based estimates incorporate the bulk of high-z star formation 
down to the faint end of the LF. We also see no evidence for 
a strong peak in the SFR versus z. This assumes that a very 
strong rise in the efficiency of producing GRBs (beyond that 
already accounted for), does not hide a drop in the SFR, al- 
though this itself would be quite interesting. While we have 
not included them in our analysis, of the three Fermi GeV- 
detected long GRBs with redshifts (shown in Fig.[TJ, two were 
at z > 3.5 (e.g., Abdo et al. 2009). Their brightness raises the 
prospect of the independent use of GeV-selected bursts. 

The current picture of small, metal-poor GRB hosts ob- 
served at low z agrees well with our GRB-inferred SFR being 
dominated by such sub-L, galaxies at high z. One might won- 
der about the whereabouts of these GRB hosts today, whether 
they continued to grow, merged into more m assive halos, 
etc. Observations of th e afterglow spectrum (e.g.. lTotani et al.l 
I2006t iMcOuinn et al.ii2008) could determine the extent that 
the host had experienced the effects of the reionizing UV 
background. Since GRBs should originate from a different 
range of overdensities than quasars, potential exists for an- 
other examination of the hierarchical history of our Universe. 

We thank Kris Stanek for helpful discussions and Nat But- 
ler for making his GRB data available and suggesting peak lu- 
minosities. We acknowledge use of the Swift public archive. 
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